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Abstract 

The current early stage in the investigation of the stability of the Kerr metric is 
characterized by the study of appropriate model problems. Particularly interesting is 
the problem of the stability of the solutions of the Klein-Gordon equation, describing 
the propagation of a scalar field of mass /i in the background of a rotating black hole. 
Rigorous results proof the stability of the reduced, by separation in the azimuth angle 
in Boyer-Lindquist coordinates, field for sufficiently large masses. Some, but not all, 
numerical investigations find instability of the reduced field for rotational parameters 
a extremely close to 1. Among others, the paper derives a model problem for the 
equation which supports the instability of the field down to a/M s» 0.97. 

1 Introduction 

The discussion of the stability of the Kerr black hole is still in its early stages. The first 
intermediate goal is the proof or disproof of its stability under "small" perturbations. 
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In the case of the Schwarzschild metric, by using the Regge-Wheeler-Zerilli-Moncrief 
(RWZM) decomposition of fields in a Schwarzschild background [26, 31, 23, 28, 24, 13], 
the question of the stability of linearized fields can be completely reduced to the question 
of the stability of the solutions of the wave equation on Schwarzschild space-time. In 
comparison to Schwarzschild space-time, the case of Kerr space-time is complicated by a 
lower dimensional symmetry group and the absence of a Killing field that is everywhere 
time-like outside the horizon. For instance, the latter is reflected in the fact that energy 
densities corresponding to the Klein-Gordon field in a Kerr gravitational field have no def- 
inite sign. This absence complicates the application of methods from operator theory, as 
used in this paper, and also of so called "energy methods" that are both employed in esti- 
mating the decay of solutions of hyperbolic partial differential equations. 

For Kerr space-time, a reduction similar to RWZM for Schwarzschild space-time is not 
known. On the other hand, the finding of a symmetry operator, containing only time 
derivatives up to first order, for a rescaled wave operator gives hope that such a reduction 
might exist [6]. If such reduction exists, there is no guarantee that the relevant equation 
is the scalar wave equation. It is quite possible that such equation contains an additional 
(even positive) potential term that, similar to the potential term introduced by a mass of 
the field, see below, could result in instability of the solutions. Second, an instability of a 
massive scalar field in a Kerr background could indicate instability of the metric against 
perturbations by matter which generically has mass. If this were the case, even a proof of 
the stability of Kerr space-time could turn out as a purely mathematical exercise with little 
relevance for general relativity. 

Currently, the main focus is the study of the stability of the solutions of the Klein-Gordon 
field on a Kerr background with the hope that the results lead to insight into the problem 
of linearized stability. Although the results of this paper also apply to the case that ji = 0, 
its main focus is the case of Klein-Gordon fields of mass /i > 0. 

Quite differently from the case of a Schwarzschild background, the results for these test 
cases suggest an asymmetry between the cases [i = and \i ^ 0. In the case of the wave 
equation, i.e., fj, = 0, rigorous mathematical, analytical and numerical results point to the 
stability of the solutions [30, 14, 10, 2, 19, 20]. On the other hand, for ^ ¥= 0, there are 
a number of analytical and numerical results pointing in the direction of instability of the 
solutions under certain conditions [11, 12, 32, 15, 18, 9, 16]. 

In particular, unstable modes were found by the numerical investigations by Furuhashi 
and Nambu for fiM ~ 1 and (a/M) = 0.98, by Strafuss and Khanna for /j,M ~ 1 and 
(a/M) = 0.9999 and by Cardoso and Yoshida for fiM < 1 and 0.98 < (a/M) < 1. The 
analytical study by Hod and Hod finds unstable modes for fiM ~ 1 with a growth rate 
which is four orders of magnitude larger than previous estimates. On the other hand, [3] 
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proves that the restrictions of the solutions of the separated, in the azimuthal coordinate, 
Klein-Gordon field (RKG) are stable for 



The paper [7], among others, gives a stronger estimate. It proves that the solutions of the 
RKG are stable for fi satisfying 



So far, these have been the only mathematically rigorous results on the stability of the 
solutions of the RKG for [i > 0. The first estimate has been confirmed numerically by 
Konoplya and Zhidenko [21]. The stronger estimate has been confirmed by S. Hod [17]. 
These results contradict the result of Zouros and Eardley, but are consistent with the other 
results above. In addition, the numerical result by Konoplya and Zhidenko finds no unsta- 
ble modes of the RKG for \iM « 1 and /iM ~ 1, contradicting above cited analytical and 
numerical investigations. 

The situation for the RKG for the case n > is quite puzzling. The RKG originates 
from the case fi = by the addition of a positive bounded potential term to the equation. 
Without the presence of a first order time derivative in the equation, from this alone, by 
help of methods from operator theory, it would be easy to prove that the stability of the 
solutions of the wave equation implies the stability of the solutions of the Klein-Gordon 
equation for non-vanishing mass. Also, the energy estimates in Lemma 4.7 of [7], indicate 
a stabilizing influence of such a term. On the other hand, so far, there is no result that 
would allow to draw such conclusion. The numerical results that indicate instability in the 
case /i ^ make quite special assumptions on the values of the rotational parameter of 
the black hole that could indicate numerical artifacts. Moreover, as mentioned before, the 
numerical investigation by Konoplya et al. does not find any unstable modes and contra- 
dicts all these investigations. Also the analytical results in this area are not accompanied 
by error estimates and therefore ultimately inconclusive. 

The present paper tries to help clearing up this situation. It derives and analyzes the 
stability of a simplified model equation from the Klein-Gordon equation, describing the 
propagation of a complex scalar field of mass /i 5= in the gravitational field of a rotating 
Kerr black hole. Also, the simplified model equation is in certain sense "mathematically 




(1.0.1) 



Here meZis the 'azimuthal separation parameter' and 



M + \/M 2 - a 2 . 
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controllable." The latter term will be given meaning in the next section. 

The remainder of the paper is organized as follows. Section 2 gives the derivation of 
the simplified model equation. Section 3 gives basic properties of operators read off from 
this approximate equation which provide the basis for the formulation of an initial-value 
problem along with a stability analysis for that equation given in the same section. The 
results of Section 3 are used in Section 4 to show the instability of the solutions of the 
approximate equation for moderately large rotation parameters a. Finally, the paper con- 
cludes with a discussion of the results and an appendix containing the proof of a result that 
is used in the main text. 



2 Derivation of the Approximate Equation 

The goal of this section is the description of the relationship between the Klein-Gordon 
equation and the approximate equation. Of course, this relationship is crucial for the model 
character of the latter equation. In particular, there are various ways to arrive at approxi- 
mate equations suitable for analytical treatment. Therefore in the following, we describe 
the derivation of the approximate equation in more detail. We mention that it is easily seen 
that our approach is in many aspects different from that of Detweiler in [12]. 

To facilitate analytical treatment, the spectral problem, that is associated with the approx- 
imate equation, is required to lead on a special case of the Coulomb wave equation. The 
latter can be reduced to an equation of confluent hypergeometric type. The solutions of 
these equations are well-known. 

We start from the Klein-Gordon equation, describing the propagation of a complex scalar 
field u e C 2 (f2, C) of mass \i 5= in the gravitational field of a rotating Kerr black hole of 
mass M > and with rotational parameter a e [0, M] in the following form 
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4 are the Boyer-Lindquist coordinates and 
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+ /i 2 E 



(2.0.2) 
= . 



A := r 2 - 2Mr + a 2 , E 



r 2 + a 2 cos 2 £ 



If not otherwise indicated, the symbols t,r,9,ip denote coordinate projections whose domains will be obvious 
from the context. In addition, we assume the composition of maps, which includes addition, multiplication 
and so forth, always to be maximally defined. For instance, the sum of two complex-valued maps is defined 
on the intersection of their domains. Finally, we use Planck units where the reduced Planck constant li, the 
speed of light in vacuum c, and the gravitational constant G, all have the numerical value 1. 
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r+ := M + V M 2 — a 2 , fl := R x (r+,oo) x (-tt.tt) x (0,tt) . 
The introduction of the new unknown u e C 2 (Q, C) by 

u := A 1/2 u (2.0.3) 
and subsequent multiplication of the resulting equation by (r 2 + a 2 ) 2 /A 2 gives 



dr 2 
1 

A 
1 d 
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1 - 
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+ /i 2 (2Mr-a 2 sin 2 6>H u 
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+4Mar 



4M 2 (2Mr - a 2 ) 



a 2 A 



(r 2 + a 2 )' 



sin 2 6» 



AA 
<5t 2 



5 2 m 



5 2 m 



(M 2 — a 2 ) > w = . 



In a first step towards the model equation, we neglect the term 



inside the factors 



(r 2 + a 2 ) 2 
a 2 A 



sin 



1 - 



sin 



(V 2 + a 2 ) 2 

multiplying the second order time derivatives of u. In addition, we neglect the term 



2 • 2 

-a sin i 



of the factor 



2Mr - a 2 sin 2 6 1 



multiplying the square of the mass of the field, /j 2 . Through this, we arrive at an equation 
whose angular derivatives and potential leave spherical harmonics invariant. The second 
step replaces the factor 

A = (r — r_)(r — r + ) , 
inside the factors multiplying terms enclosed by curly brackets, by 

{R-rJ){r - r+) , 
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where R e [r + , oo) is a new parameter. Finally, the third step replaces the r inside terms 
enclosed by curly brackets by the constant map of value R. Through the last two steps, we 
arrive, after separation of the t, 9, and ip variables, at a special case of the Coulomb wave 
equation. The latter can be reduced to an equation of confluent hypergeometric type. The 
solutions of these equations are well-known. Through the previous steps, we arrive at a 
form of the model equation 



= - 



+ 



+ 



dt 2 



(2.0.4) 
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AM 2 (2MR - a 2 ) — + AMaR 



ct 2 



sm 2 9 d(p 2 
d 2 u 



+ 2MRf/ 



dtdp dip 



- (M - a 2 ) 



R- 



R-r- 



We note that i?_ = if and only if R = r + and a = M. For this reason, in the following, 
we implicitly assume for the case R = r + that a e [0, M). Substitution of (2.0.3) into the 
(2.0.4) leads to the model equation for (2.0.2) 



1 + 



AM(R + M) 4M 2 (2MR-a 2 ) 



+ 



R 2 _{r-r+f 



+ 



+ - 
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We note that 



+ (M 2 - a 2 ) 
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M{R + M) > M{r+ + M) 



2MR - a 2 > 2Mr+ - a 2 



R 2 _ 
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66 
- 1 



sin 2 9 dip 2 



+ 2MRff 



+ ix 1 \ u = 



(2.0.5) 



(r 2 + a 2 + 2M 2 ) > , 



rt > 



Division of (2.0.5) by 



p ._ 1+ 4M(R + M) + AM 2 (2M R - a 2 ) 



R-(r- 



R 2 _ (r-r_, 



R 2 _{r - r+) 2 + AM(R + M)i?_(r - r+) + 4M 2 (2MR - a 2 ) 
R 2 _{r-r+) 2 

[R_{r - r+) + 2M(R + M)f - AM 2 (R 2 + M 2 + a 2 ) 

R 2 _{r-r+) 2 ' !1 
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gives 

d 2 u 
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AMaR 
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R±(r-r+) 2 
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+ [M 2 - a 2 ) 



R 2 



(r — r_) 2 



+ p 2 \ u = 



(2.0.6) 



Substitution of the ansatz 



u(t, r, 9, ip) = v(t, r)Yi m (6, ip) 



where v e C 2 (R x I,C), I := (r+,oo), to e Z and Z e {|m|, |m| + 1, . . .}, c; := Z(/ + 1) 
and Y/ m the spherical harmonics corresponding to m and Z, into (2.0.6) gives the reduced 
equation 



d 2 v , 8v 



Here 



4maM R 
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pA 5r 



4maMR 



(2.0.7) 



pi^O - r + ) 2 i?2(r - r + ) 2 + 4M(i? + M)R_{r - r+) + 4M 2 (2MR - a 2 ) 



V 



o c t + 2MRu 2 
p + 



{M 2 -a 2 )[l-^] 



+ m 2 a 2 



i?_(Y-r + ) i?2_(r-r + ) 2 
M 2 i£(r - r + ) 2 + (q + 2MRp 2 )R^(r - r+) - [{M 2 - a 2 ) [l - j—p] + ra 2 a 2 ] 



R-{r - r+) 2 + 4M(i? + M)i?_(r - r+) + AM 2 (2MR - a 2 ) 



We note in particular that b and 1/ are bounded continuous real-valued functions that have 
continuous extensions to /. We define 



V : = inf V 

re I 



and note that 



m 2 a + M 2 — a 2 
AM 2 {2MR- a 2 ) 
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Fig. 1: Sketch of Graph(V) from r + to 10, for M = 1, a = 1/2, m = 3, R = r+, A* = 0.5 and Z = 20. 



3 Formulation of an Initial Value Problem for the Re- 
duced Equation 

This section gives basic properties of operators read off from the approximate equation 
(2.0.7) as well as the formulation of an initial-value problem along with a stability analy- 
sis for the latter, following [4, 5]. 

As data space for (2.0.7), we choose the weighted L 2 -space 

X:=L&I,pA) . 

Since 

. 4M(R + M), , AM 2 (2MR- a 2 )(r -r\ 

p A = A + K — ; (r - r_) + K — 2 ^ '- 



R±(r-r+) 



this implies, in particular, that all continuous functions that have continuous extensions to 
/ which are non-vanishing on r + are not part of the data space. 

Definition 3.1. We define B as maximal multiplication operator in X corresponding to b. 

Lemma 3.2. B is a bounded linear and self-adjoint operator on X. B is positive for 
rn^O and negative for m 0. 

Proof. Since b is in particular bounded an real-valued, B is a bounded linear and self- 
adjoint operator on X. Further, since b is positive, negative for m 5= and m ^ 0, 
respectively, B is positive, negative for m > and m 0, respectively. □ 
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Definition 3.3. We define the linear operator Aq in X by 

1 



A .f:=-—(Af')' + Vf 



for every / e Cq (I, 1 



Lemma 3.4. A is a densely-defined, linear, symmetric, semibounded and essentially 
self-adjoint operator in X with lower bound Vq. 

Proof. First, we note that by 

Uf := (pA)" 1 / 2 / 
for every / e £ 2 (/), there is defined a Hilbert space isomorphism 

U:L&I)^Ll(I,pA) . 

In particular, the operator 

A ou := U-'AqU 
in Lg(I) has the domain Cq (I, M) and is given by 

A 0C // = -(/'/p)' + (^ + Vc / )/ 

for every / 6 Cq (I, M), where 

Vu := -- (A' 2 + 2AA")F- 1 + - (2A'F' + AF")F- 2 - - A 2 F' 2 F- 3 , 
4 2 4 

F := P A 2 . 

We note that Vu is in particular bounded. Hence it follows also that V + Vu is bounded. 
We define the auxiliary operator A 0U t m ^c(-0 by 

AouTf:=-(f'/p)' 

for every / e Cq (J, C). Since Cq (I, C) is dense in £ 2 (/), Aqut and A t/ are densely- 
defined. Also, Aqut and ^4qc7 are obviously linear. Further by use of integration by parts 
along with the boundedness of V + Vu, it follows that Aqut, Aqu are symmetric and 
semibounded from below. The equation 

-(/7p)' = o 

admits the C 2 -solution 

4M(fl + M) 1 . , 4M 2 (2MR-a 2 ) 
/==r-r + + — ^-^ln(r-r + ) ^ (r _ r+) J - 
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Since / is not in L^(I) at both ends of I, it follows that Aqjjt is in the limit point case 
at r + and at +00. Hence Aqjjt is essentially self-adjoint (see, e.g., [29]). Further, since 
V + Vjj is bounded and real-valued, it follows from the Rellich-Kato theorem, e.g., see 
Theorem X.12 in Volume II of [25], that Aqjj is also essentially self-adjoint and that the 
domains of the closures Aqjjt of Aqjjt and Aqjj of Aqjj coincide. Hence, we conclude 
that Aq is a densely-defined, linear, symmetric, semibounded and essentially self-adjoint 
operator in X. Finally, it follows from integration by parts that Vq is a lower bound for 



Definition 3.5. In the following, we set A := Aq. 

From Lemma 3.4, it follows that 

Theorem 3.6. A is a densely-defined, linear, symmetric, semibounded and self-adjoint 
operator in X with lower bound Vq. 

Further, from the proof of Lemma 3.4, we can conclude the following, see, e.g., [29]. 
Corollary 3.7. 



Aq. 



□ 



(i) For every z e C\(M x {0}), / 6 L|(7, pA) and re I, 

[{A-z)- 1 f]{r) = -v 2 {r) f Vl (r')(pAf)(r ')dr' 





for some cel. 



(ii) By 



Aq*v := 



1 



(Av 1 ) 1 + Vv 



pA 
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As a consequence of the self-adjointness of A and semiboundedness of A with lower 
bound Vq, the objects X, A-v B +s ■= A—Vo + £,B and C : — (— Vq + e) are easily seen 
to satisfy Assumptions 1 and 4 of [4]. 1 Here e > is assumed to have the dimension l~ 2 . 
The exact value of e does not affect subsequent results in any essential way. Application 
of the results of [4] give, in particular, the following well-posed formulation of the initial 
value problem for (2.0.7). 

Theorem 3.8. (Formulation of an initial value problem for the reduced equation) 

(i) By 

Y := D(A% +e ) x X 

and 

m := {A% j+E b\A% +e m ) + <6l%> 
for all £ = , £,2) , V = (^lj ^72) e Y, there is defined a complex Hilbert space 

(y,(\))- 

(ii) The operators G and — G defined by 

G(d,v) ■= (-mAZ + iBr,) 
1/2 

for all £ e D(A) and 77 e D(A_ V +e ) are infinitesimal generators of strongly contin- 
uous semigroups T + : [0, 00) — ► F) and T_ : [0, 00) — ► L(Y, Y), respectively. 

1/2 

(iii) For every to e M and every £ e -D(A) x D(A_ V +£ ), there is a uniquely determined 
differentiable map u : K — » 1" such that 

u(to) = £ 

and 

u '(t) = -Gu(£) (3.0.8) 

for all t e 1, Here ' denotes differentiation of functions assuming values in Y. 
Moreover, this u is given by 

1 j • { T_(-£)£ fori <0 1 J 

for alii eK. 

(iv) For all £ e [0, 00): 

|r ± (£)Kexp(||C||||A:^ 2 +£ ||£) , 
where | |, || || denote the operator norm for L(Y, Y) and L(X, X), respectively. 

1 See also the Section 5.1 on 'Damped wave equations' in [5]. 
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(v) A e C is a spectral value, eigenvalue of iG if and only if 

A — XB — X 2 
is not bijective and not injective, respectively. 

(vi) For any A from the resolvent set of iG and any 77 = (771, 772) 6 Y one has: 

(iG-X)- 1 r 1 = (fci^ + m)) , 

where 

e := (A - XB - A 2 )- 1 [(B + X) m - i m ] . 

Equation (3.0.8) is the interpretation of (2.0.7) used in this paper. In this sense, (iv) shows 
the well-posedness of the initial value problem for (2.0.7), i.e., the existence and unique- 
ness of the solution and its continuous dependence on the initial data. Moreover (3.0.9) 
gives a representation of the solution and (iii) gives a rough bound for its growth in time. 
In general, this bound is not strong enough to show, if applicable, stability of the solutions 
to (2.0.7). Part (v) reduces the determination of the generator's spectrum to the determina- 
tion of the spectrum of the operator polynomial A — XB — X 2 , A e C [22, 27]. Moreover 
(vi) does the same for the resolvents. 

Lemma 3.9. For every A e C, 

ker(A — XB - A 2 ) = ker(A * - XB - A 2 ) . 
Proof. Let A e C and / e ker(A — XB — A 2 ). Then 

f = (A + i)- 1 (X 2 + XB + t)f. 
Hence the statement follows from Corollary 3.7 and the fact that b is a C^-function. □ 

Assumption 3.10. In the following, we use the conventions that yfz, for some z e C, 
denotes some square root of z, whereas [ j 1 / 2 : C\(— 00, 0] — > C denotes the principal part 
of the complex square root function, 

Theorem 3.11. If A e C\{— fx, 7 e C satisfies 

7 2 - 7 + [4M 2 (2Afi? - a 2 ) A 2 + AmaMRX + m 2 a 2 + M 2 - a 2 ] = , (3.0.10) 
then the C 2 -solutions of 

- ^-(Av')' + (V - Xb- X 2 )v = (3.0.11) 
pA 
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are given by 

A" 1 / 2 (r - r+y< e vV-* 2 r w (-2^fi 2 - A 2 (r - r+)) , 
where w : £1 — * C is some solution of Kummer's equation in the complex domain 

zw"{z) + (27 — z)w'(z) — aw{z) , 
z e U, where U c C is an open neighborhood of the half-ray 



— V yU 2 — A 2 . (0, go) 

and 

AM(R + M)\ 2 -ci- 2MRu 2 
a := 7 H . . 

Proof. Direct calculation. □ 

Hence, it follows from known properties of the solutions of the confluent hypergeometric 
equation that 

Corollary 3.12. If A e C\{-(i,[i}, 7 e C satisfies (3.0.10), and 2 7 $ Z, then all C 2 - 
solutions of (3.0. 1 1) are linear combinations of the linear independent solutions 

A- 1 / 2 (r - r+ y eV " 2 ~ x2 r M{a, 2 7 , -2vV 2 - A 2 (r - r+)) , 

A" 1 / 2 (r - r+) 1 -T e vV-A 3r W( 1 + a - 2 7 , 2(1 - 7), -2yJ y? - A 2 (r - r+)) , 

where for /3 e C\Z, M(a, /3, •) denotes the confluent hypergeometric function of the first 
kind, as defined in [1]. 

For A e C, (3.0.10) has the solutions 



11/ R 2 
7 = - + — \j-AM 2 {2MR - a 2 )\ 2 - AmaMRX - m 2 a 2 - (M 2 - a 2 ) + — 



(3.0.12) 



7 = - + n/ ^ 2 _ a2 , 1/2 (2Mr+X + ma) . (3.0.13) 



In the special case R = r+, the latter reduces to 

1 i 

2 ~ 2(M 



Hence in this special case, the condition 27 ^ Z is equivalent to the condition 

A ^ 1 — [ma + ik(M 2 - a 2 ) 1/2 l , (3.0.14) 

2Mr + 

13 



00 (c 



where k e Z. Further, straightforward calculation shows that argument of the square root 
in (3.0. 12) is in (— oo, 0] if and only if A e S, where 

-maR- \m 2 a 2 A(R) + (2MR-a 2 )\M 2 -a 2 - (R 2 _/A)\] 1/2 

— 00 

2M(2MR- a 2 ) 

maR+[m 2 a 2 A(R) + (2MR - a 2 )\M 2 - a 2 - (i? 2 /4)|] 1/2 
2M(2MR - a 2 ) 

In particular, we note for the special case R = r + that 

S = R . 

Hence 

Lemma 3.13. For A e C, (3.0.10) has a solution with real part > 1/2 if and only if 
A e C\S. In the remaining cases, this equation has only solutions with real part 1/2. If 
A e C\S, the solution with real part > 1/2 is given by 



1 1 

2 + Rl 



-AM 2 (2MR - a 2 )X 2 - AmaMRX - m 2 a 2 - (M 2 - a 2 ) + 



R 2 ^ 2 



(3.0.15) 



and the other solution of (3.0.10) has real part < 1/2. 
Hence, we conclude from Corollary 3.12 that 

Corollary 3.14. If A e C\{— fx, fi}, and the solutions of (3.0.10) are not part of (1/2).Z, 
then there is a non-trivial solution v to (3.0.1 1) such that x (r c) ■ v e X for some c > r + 
if and only if A 6 C\S. In the latter case, all such v are multiples of 



A" 1 / 2 (r - r+y e V>^ r M(a, 2 7 , -2^J - A 2 (r - r+)) , 
where 7 is given by (3.0.15). 

If A e C\((-oo, -/x] u \p, 00)), and the solutions of (3.0.10) are not part of (1/2) .Z, it 
follows from the asymptotic properties of the confluent hypergeometric function of the 
first kind, see, e.g., [1] 13.1.4, that 

A-V2 ( r - r+ y e-^ 2 -^ 1 ' 2 r M(a, 2 7 , 2(^ 2 - A 2 ) 1 ^ ( r _ r+ )) , 

is exponentially growing for r — * 00 if a $ — N. Hence, by use of Lemma 3.9, we arrive at 

Theorem 3.15. If A e C\((— 00, — ^] u [fi, 00)), and the solutions of (3.0.10) are not part 
of (1/2).Z, then ker(A - XB - A 2 ) is non-trivial if and only if A e C\S and 



1 1 

2 + RZ 



R 2 

-4M 2 (2MR - a 2 )X 2 - AmaMRX - m 2 a 2 - (M 2 - a 2 ) + — 



1/2 
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for some n e N. In the latter case, ker(A — XB — A 2 ) consists of the multiples of 

A" 1 / 2 (r - r + y e-^ 2 -^ 1 ' 2 r M(-n, 2 7 , 2( M 2 - A 2 ) 1 / 2 (r - r+)) , (3.0.17) 
where 7 is given by (3.0.15). 

Since in the case that /i 2 [2M(R+ 2M )] _1 Q, the left hand side of (3.0.16) has real part 
> 1/2, whereas the right hand side has real part 0, we obtain the following corollary. 

Corollary 3.16. If 

2 1(1+1) 



2M(R + 2M) 

A e C\M, and the solutions of (3.0.10) are not part of (1/2) .Z, then kev(A - XB - A 2 ) is 
trivial. 

In particular, in the case of vanishing mass fj, = and for the considered cases, this implies 
that kcr(A — XB — A 2 ) is trivial. Also, we obtain the following estimate for the case that 
the previous inequality is violated. 

Corollary 3.17. If 

, 1(1 + 1) 



M > 



2M(R + 2M) 

X e C\M, the solutions of (3.0.10) are not part of (1/2).Z and kei(A - XB - A 2 ) is 
non-trivial, then 

|A| < + ' 2M ' R + 2A y-'" + "1 2 }' /2 ■ (3.0.!8) 

Proof. Straightforward estimates yield that the inequality 

a a [2M(R + 2M)^-l(l+l)f 
l A ~ A 1 I > ^5 

implies that 

i \2M(R+2M)ii 2 - q] ^ () 



- - Re 
2 



2i?_(^ 2 -A 2 )V2 

and hence that the left hand side of (3.0. 16) has real part > 0, whereas the right hand side 
has real part ^ 0, which contradicts that ker(^4 — XB — A 2 ) is non-trivial.^ Hence, 

a 2 | ^ [2A/(fl+2Af)/x 2 -Z(/ + l)] 2 
|A ~ M|< i£ ' 

The latter inequality implies (3.0. 18). □ 
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Lemma 3.18. If A is such that the solutions of (3.0.10) are not part of (1/2). Z, then 
A e K x (—oo, 0) is such that ker(A — XB — A 2 ) is non-trivial if and only if 



A 



2z 



for z e (C\B M (0)) n ((0, go) x R) satisfying 



M\M A + a 2 + R z )z* + 2M 



n+ g ) R-(R+ M) nnaf! 



+ [(M 2 - a 2 ) + m 2 a 2 + n(n + 1)R 2 _ + 2M(R + M)(cj - 2M 2 /x 2 )] 



+ 2 •{ ( n + - ) iL. [ci + M{R - M)n 2 ] - imaMRfi 



+ | cf + 2 [(M 2 - a 2 ) + m 2 a 2 + n(n + 1)^ - 2M 2 q] /i 2 



+ 2M 2 {R 2 + AMR + 3M 2 - a 2 )/i 4 jz 4 

n + ij -R- [ci + M(R— M)fi 2 ] + imaMR/j, 2 \ [i 

,2„2 



2 z 3 



+ [(M 2 - a 2 ) + m 2 a 2 + n(n + 1)R 2 _ + 2M(R + M)(q - 2/i 2 M 2 )] /x 4 z 2 
1 



+ 2M 
and such that 



n + - ) R-(R + M) + imaR 



^z + M\M Z + a z + R 2 )fi a = 



1 



2 2 

n + ^ I R- - M(R + M) Z +M - [q - 2M(R + 2iW> 2 ] „ Z . 

2 / Z ' Z z + [L A 



is part of (0, oo) x E. 

Proof. The statement follows by direct calculation from Theorem 3.15 and Lemma 6.1. 

□ 

Lemma 3.19. If R = r+, i.e., i?_ = 2(M 2 - a 2 ) 1 / 2 , and A satisfies (3.0.14), then 
A e K x (—oo, 0) is such that kcr(yl — XB — A 2 ) is non-trivial if and only if 



A 



for z e (CYB^(O)) n ((0, oo) x R) satisfies 



4 (2n+ 1)(M 2 -a 2 ) 1 / 2 + ima 3 q - 2M 2 /i 2 2 
M(2r+ + M) - z + M ( 2r+ + M) Z 
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(2n + l)^-a^ + ima M_ 
M{2r+ + M) p 2r+ + M p 

Proof. The statement follows by direct calculation from Theorem 3.15 and Lemma 6.1. 

□ 

We note that if z £ (C\B^(0)) n ((0, oo) x R) satisfies (3.0.19), then z* is a solution of 
(3.0.19), where m is replaced by — m, that is contained in (C\i? M (0)) n ((0, oo) x R). 



4 Unstable modes 

In this section we proceed to solve (3.0.19) and show the existence of unstable modes in 
a subregion of the parameter space. Although the roots of any forth degree polynomial, 
such as that in (3.0.19), are known explicitly, the expressions, seen as functions of all the 
parameters involved (i.e. a, fi, n, I and m) 1 , are too complicated to get any intuitive under- 
standing of the problem just from their analytical form. Thus, we find it more insightful 
to simply do a direct numerical evaluation at different parameters' values. This approach 
allows us to find the unstable modes described below. 



Before continuing, we mention a few words about the numerical errors. Besides con- 
ducting error estimations of the eigenvalues found, we also checked our results by finding 
the roots with the Newton Raphson (complex) method and then compare them with the 
values found by direct evaluation, giving identical results up to error. We find that the 
errors can be made as small as desired just by increasing the number of figures used in the 
calculations. All values presented in this paper are accurate at least up to the last figure 
shown. 



We find unstable modes in a range of values of fi, which increases with I, m and a. That 
region decreases with n, but the dependence on n seems to be much smaller than that on 
I, m and a, and hence for the majority of the results presented here we just set n = 0. 
Let the reader assume n = from now on unless otherwise specified. We find unstable 
modes for I and m starting at 3 and large a, mostly with I = m, although for large I we 
also find unstable modes for some m < I. To give a few examples, for I = m = 3 we 
find instabilities with a = 0.9976 and larger, while for I = m = 10 we find them for 
a = 0.9833 and larger. A more complete description is given in the figures. The curves 
in Figure 2 represent unstable regions in [i corresponding to four values of a, for fixed 
I = in; and to four values of / = m and fixed a. In each case, we show the complete range 
of values fi for which we found unstable modes (vertical axis of the figures). Interestingly, 
in all the cases studied (including others not shown in the figures) that range starts exactly 
(up to numerical accuracy) at the limit given in Corollary 3.16. Finally, in Figure 3 we 

1 Since the solutions are scalable in M we will only present results for M = 1. 
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Fig. 2: Eigenvalues A of unstable modes (A < 0) for varying p, a and I = m. /J, is shown in the vertical 
axis, while the horizontal axes represent the real and imaginary part of A. The dashed lines are a projection 
of the curves into the complex plane to improve visibility, left panel: I = m = 10; a = 0.984 (red), 0.986 
(green), 0.988 (blue) and 0.99 (magenta), right panel: a = 0.99; I = m = 6 (red), 8 (green), 10 (blue) and 12 
(magenta). 

show the norm of the eigenfunctions (3.0.17) (see Theorem 3.1 1) as a function of r — R + 
for some values of [i. In the inset we use a logarithmic scale in the horizontal axis to show 
the details close to R + . 

Since previous works have reported instabilities only for large values of a, it is worth 
searching for instabilities with the smallest possible a. Figure 4 shows the minimum value 
of a, a m i n , for which we could find instabilities for given I, For this task we only investi- 
gated modes with m = I, since usually setting m < I gives rise to less (if any) unstable 
modes. Although an exhaustive search was made to determine a m j n , it is impossible to 
strictly rule out the possible existence of unstable modes beyond a m j n with our method 
of direct evaluation, since an infinite number of such evaluations would be necessary to 
completely cover any subregion of the parameter space (unless the continuous parameters, 
a and /i, are both fixed). In general, for given a a range of /i values give unstable modes. 
If a is then decreased, that range gets smaller until, at a = a mm , only one value of ^ gives 
unstable modes (see the left panel of figure 2). Each of these values are show as annota- 
tions in Figure 4. As seen in the figure, a min decreases with / and seems to asymptote to a 
value a asym m 0.9718. It would be interesting to determine the mathematical root of this 
threshold. This task will be pursued in a future work. 
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Fig. 3: Eigenfunctions of unstable modes with n = 0, 1 = m = 10, a = 0.988; and /i = 4.184 (solid), 4.237 
(dash), 4.290 (dot dash), 4.396 (long dash), 4.502 (dots). 



0.995 - 



= 0.985 - 



0.975 - 



0.97 V 



t V 195 



.3 



.31* 
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100 
1, m 



1000 10000 



Fig. 4: Minimum a for which unstable modes were found for given I = m. The annotated numbers show the 
approximate value of ^. 
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5 Discussion of the Results 



From the Klein-Gordon equation in a Kerr background, we derive a simplified model equa- 
tion showing the instability of the reduced, by separation in the azimuth angle in Boyer- 
Lindquist coordinates, field for a certain range of the masses [i of the field. In addition, 
we give a well-posed initial value formulation of that reduced equation, along with a sta- 
bility analysis of the corresponding field. The latter shows instability down to rotational 
parameters a/M <%j 0.97. This result supports claims of previous analytical and numeri- 
cal investigations that show instability of the massive Klein-Gordon field for a extremely 
close to 1. 

From here, mathematical investigation could proceed in 2 directions. First, it might be 
possible, to use the model for the proof of the instability of the massive Klein-Gordon 
equation in a Kerr background, using a perturbative approach. Another direction consists 
in further analysis of the model in order to find the mathematical root of the instability as 
well as an abstraction to a larger class of equations that includes the massive Klein-Gordon 
equation on a Kerr background. 
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6 Appendix 

Lemma 6.1. Let (i e [0, co). By 

for every A e E x (— oo, 0), where (fi 2 — A 2 ) 1 / 2 denotes the square root with strictly 
positive real part, there is defined a biholomorphic map 

g :R x (-oo,0) — (C\S M (0)) n ((0,oo) x E) 

with inverse 

g' 1 : (C\B M (0)) n ((0, oo) x E) — E x (-oo, 0) 
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Im(z) 




Fig. 5: Sketch of domain and range of the map g from Lemma 6. 1 . 

given by 

9-\z) = - l(z* - M 2 ) 
for every z e (C\B M (0)) n ((0, oo) x E). In addition, 

for every z e (C\B M (0)) n ((0, oo) x E). 

Proof. First, we note that g is well-defined. For this, let A e E x (— oo, 0) and A 
Re(A), A 2 := Im(A) < 0. Then 

/i 2 - A 2 = fi 2 - \\ + \ 2 2 - 2UiA 2 . 

Hence /i 2 — A 2 is real iff Ai = 0. In the latter case, 

Re (/i 2 - A 2 ) = n 2 + A 2 . > . 

Therefore, /i 2 — A 2 e C\((— oo, 0] x {0}), and there is precisely one square root of fi 2 
with strictly positive real part. Further, 

Re [iX + (jit 2 - A 2 ) 1 / 2 ] = -A 2 + Re ( M 2 - A 2 ) 1 / 2 > . 

In particular, if 

z := iX + (/i 2 - A 2 ) 1 / 2 
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is such that \z\ < fi, then 

z 2 - 2iz\ - A 2 = (z - ^A) 2 = M 2 - A 2 . 

Hence 



z 2 — = 2izX 



and 



The latter implies that 



2 V z 



Xn = ImA = 



Re(z) 



and hence that Re(z) O.i As a consequence, |z| > /x. 



For the second step, let z e (C\B P (0)) n ((0, oo) x E), a; := Re(z), and j/ := Im(z). 
Then 



-2 [Z 



fir 

z 



2/1 + 



and 



— - I a; + ty — 



— ix I 1 — 



x + ly 



e E x (-oo,0) 



+ 



^ 4 I z 



z — 



z + 







z / 


2\ 2 


1/2 






T 





1/2 



Since 



_ 1 
~ 2 

the latter implies that 



a; + ij/ + -j— (x - iy) 



x 1 + fio )+ *V 1 



e (0,co) x 



5 -o *~ 



1/2 
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In particular, we conclude that by 



h(z) :- 



2z 



for every z e (C\B M (0)) n ((0, oo) x R), there is defined a map 

h : (C\B M (0)) n ((0,oo) x E) — E x (-00, 0) 

such that 

(goh)(z) = z 

for every z e (C\-B M (0)) n ((0, 00) x E). Therefore, 5 is surjective. 



Further, if /x = 0, 



(h o 5 ) (A) = - - [iX + (-A 2 ) 1 / 2 ] = _ _ (iA + iX ) = A 



and if fi ^ 0, 



(fco 5 )(A) 



iA+G^-A 2 ) 1 / 2 - 



iA + Gu 2 -A 2 ) 1/2 -M 2 



zA + (it 2 - A 2 )V2 
iA-(/i 2 -A 2 )V2 



A 2 - (it 2 - A 2 ) 
*A+( M 2 -A 2 ) 1 / 2 + zA-( M 2 -A 2 ) 1 ^ 



A 



for every A e E x (-co, 0). Therefore, g is injective and altogether bijective with inverse 
h. Finally, it follows for A e E x (-co, 0) that 



- ( M 2 - A 2 ) 1 / 2 = iX - g( X ) = ih(g(\)) g(X) = -i l - (g(X) 



9(A) 



-ff(A) 



gWJ 2 5 (A) 



[(9(A)) 2 + M 2 ] 



and hence for every z e (C\B M (0)) n ((0, 00) x E) that 

[mM<tW] 1/2 = ^(* 2 + m 2 ). 



□ 
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